ASEP on a ring conditioned on enhanced flux 
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We show that in the asymmetric simple exclusion process (ASEP) on a ring, 

conditioned on carrying a large flux, the particle experience an effective long- 
> , 

range potential which in the limit of very large flux takes the simple form U 
ON 
OO 

r**"- ' ilar to the effective potential between the eigenvalues of the circular unitary ensemble 

in random matrices. Effective hopping rates and various quasistationary probabil- 
ities under such a conditioning are found analytically using the Bethe ansatz and 

• i— 1 . 

Kn ' determinantal free fermion techniques. Our asymptotic results extend to the limit of 

5-j ' 



-2 Yli^j l°g I smnirii/ L — rij / L)\, where m, n<i, . . . n^ are the particle positions, sim- 



large current and large activity for a family of reaction-diffusion processes with on- 
site exclusion between particles. We point out an intriguing generic relation between 
classical stationary probability distributions for conditioned dynamics and quantum 
ground state wave functions, in particular, in the case of exclusion processes, for free 
fermions. 

PACS numbers: 05.40.-a, 02.50.Ey, 82.20.-w 



I. INTRODUCTION 

In this work we investigate the spatio-temporal structure of a family of stochastic inter- 
acting particle systems, specifically exclusion processes, in the case of the extreme event in 
which the system supports a very large current for an extended period of time, i.e., while 
the system is in a state very far from typical behaviour. As this type of non-equilibrium 
behaviour does not have a long history of study in statistical mechanics, we first make some 
brief remarks on how this problem is situated in general terms in the study of the statisti- 
cal properties of extreme events. Then we state the precise problem that we are going to 
address. 

When studying statistical properties of extreme events [l[, i.e., of highly non- typical 
values of some physical observable in a many-body system with noisy dynamics, the first 
question to ask is the full probability distribution of that observable since one wishes to 
know not only fluctuations around its mean, but also the distribution very far from it. 
From that knowledge one may be able to infer knowledge on the extreme value statistics 
for that observable. Then, a second problem of interest is the spatio-temporal patterns that 
generate such atypical values of observables over some period of time. Even though very 
difficult to tackle, this is an important question since a good understanding of this aspect 
of non-equilibrium behaviour, which lies beyond a purely statistical description in terms of 
extreme value statistics, may give insight into the physical mechanisms that generate extreme 
events. One could potentially identify and apply forces to enhance the probability of such 
rare events. Notice that in this setting we are not interested in searching for just any force 
that would make some rare large fluctuation of a physical observable typical but for those 
very specific forces that retain the spatio-temporal patterns of the original unforced process 
that generates these rare fluctuations by its own intrinsic random dynamics. Conversely, if 
suitable forces would suppress the relevant space-time patterns, one could make such events 
even more unlikely. 

Carrying out such an analysis is obviously very difficult even for relatively simple complex 
systems. However, this problem turns out to be interesting already in the context of simple 
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exactly solvable statistical mechanics models. To be specific, we consider in this work the 
Asymmetric Simple Exclusion Process (ASEP) [2|, y| and we look for the spatio-temporal 
patterns that generate an atypically large flux over a long period of time and for which 
effective hopping rates (arising from some extrinsic applied force) would generate such a 
state as typical stationary state. Earlier exact work on currents below the typical value 
have revealed that long-lasting low-current events are supported by travelling waves which 
represent density profiles with two regions of constant densities o-i o | l2|. S uch spontaneous 
phase separation behaviour has earlier been suggested by Fogedby [13|, [14J who used renor- 
malization group arguments to show that atypical low-current fluctuations are realized by a 
gas of shock and antishock solutions of the noisy Burgers equation that describes the large 
scale behaviour of the ASEP. In the context of phase separation these shocks and antishocks 
are the domain walls that separate regions of high and low density. On a microscopic level 
the shocks are known to be sharp, the bulk properties of each phase are reached quickly 
within a small distance from the shock position |4j. 

Understanding how the particle system organizes itself microscopically for atypically large 
currents, and how such a microscopic structure could be generated as typical event by intro- 
ducing suitable forces, requires a different and careful exact analysis. A step in this direction 



was done in recent work by Simon [15], where the effective hopping rates of the particles 
in specific regimes of the ASEP with open boundaries were computed. To make further 
progress and to investigate an arbitrary number of particles we consider in the present work 
periodic boundary conditions. Using Bethe ansatz we compute the quasi-stationary dis- 
tribution of the system conditioned on sustaining a large current and we obtain an exact 
effective interaction potential between particles that would make this quasi-stationary dis- 
tribution stationary, i.e. typical. We also stress an interesting analogy with random matrix 
ensembles. 

It will transpire that for very large currents our results for the ASEP are rather generic. 
The results remain valid for all exclusion processes with nearest neighbour hopping, irrespec- 
tive of the hopping asymmetry or even particle number conservation. Hence our approach 
covers also e.g. react ion- diffusion processes of the type studied in 5|, |6j. These include the 
well-known contact process, diffusion-limited annihilation, the coagulation/decoagulation 
model and other well-studied processes. Moreover, the same results are valid for the limit of 
large hopping activity, irrespective of direction. We also point out an intriguing connection 



to quantum mechanical expectation values for free fermion systems. 

The paper is organized as follows: In the following section (JTTJ) we define the problem and 
present our approach. Mainly for pedagogical reasons this is briefly illustrated for the case 
of just two particles (Sec. IIIIj) and then in detail for three particles in Sec. IIVI From these 
preliminaries we proceed to the investigation of an arbitrary number of particles in the limit 
of very large current in Sec. |V] where our main results are derived. We finish with a brief 
summary and discuss generalizations of our results to other particle systems (Sec. IVIj) . 



II. ASEP CONDITIONED ON CARRYING AN ATYPICAL CURRENT 

In the asymmetric simple exclusion process (ASEP) [2|, |3[ in one dimension with peri- 
odic boundary conditions and L sites, particles jump independently after an exponentially 
distributed random time with mean l/(p + q) to a nearest neighbor site, provided that the 
target site is empty. The probability to choose the right neighbour (clockwise in a periodic 
chain) is p/(p + q) while the probability of choosing the left neighbour (anticlockwise) is 
q/(p + q). Throughout the paper, we assume p + q = 1 for notational simplicity. If the 
target site is occupied the jump attempt is rejected. This exclusion principle guarantees 
that each site is always occupied by at most one particle. An instantaneous configuration 
of this system with N particles can therefore be represented by an ordered set of integer 
coordinates {n} = {ni, . . . , rijv} where rik+i > n\. and n^ G {1, . . . , L} for all k. 

Mathematically, the dynamics of the ASEP can be defined through a master equation for 
the probability P m (n,t) to find a configuration n at time t starting from a configuration m 
at time 0. The master equation takes the form 

— P m (n, t)=^2 [w ntU >P m (ri, t) - uv,„P m (n, *)] , (1) 

where "UV,n is the transition rate (0, p or q for the ASEP) to go from a configuration n 
to a configuration n . By integrating ([1]) one obtains the solution of the master equation 
for any given initial configuration m, i.e., the conditional probability to find a particle 
configuration n at time t, given that the process started from configuration m. In the 
stationary distribution of a system with N particles, all particle configurations are equally 
likely. Therefore, up to finite-size corrections of order 1/L, the stationary particle current 
has the simple form j = (p — q)p(l — p) for particle density p = N/L. For details, see 



e.g. UM\. 



Since the master equation can be used generally for describing Markov processes 
we shall use at some places below the generic symbol C for microscopic configuration of a 
process. 

From the linear form of ([1]) it is obvious that the master equation can be cast in matrix 
form 

j t \P{t)) = -H\P(t)} (2) 

where the probability vector \P(t)) has all the probabilities P(n,t) as its canonical compo- 
nents. In a conveniently chosen tensor basis the stochastic generator H of the ASEP, often 
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called "quantum Hamiltonian" , takes the form [3] 

L 

H = -H [P(4 s k+i - Ml - nk+i)) + qfa s+ +1 - (1 - n k )n k+1 )] = W + + W. + W (3) 
fc=i 

Here s k = (<t^±<t^)/2 are the SU(2) spin-1/2 ladder operators acting on site k of the lattice 
and rik = (1 — <Tfc)/2 is the projection operator on states with a particle on site k. The 
operators W + , W-, and W correspond respectively to the pure right-jump generator, to the 
pure left-jump generator and the diagonal part (see below for a detailed discussion). The 
stationary distribution is encoded in the right eigenvector of H corresponding to the lowest 
eigenvalue 0. 

The process as defined above describes only the evolution of the microscopic particle 
configuration, but does not keep track of the current that flows during the evolution. This 
can be achieved by introducing another random variable which counts the number of jumps 
across a given bond k, k + 1 of the lattice, see [7j for a detailed description. Then each 
time a particle jumps across that bond to the right (left) the value of the current counter 
is incremented (decremented) by one unit. After some time t this random variable then 
provides the integrated current, i.e., the total net number of jumps </&(£) across bond site 
(k, k + 1) up to time t. Because of particle number conservation one may equivalently one 
may also count the total net number of jumps in the lattice. We denote this quantity by 
J(t). By the law of large number one expects an asymptotically linear growth J(t) oc jt 
where j is the stationary current of the ASEP. 

With these definitions at hand we can sharpen the question of space-time realizations 
for extreme currents and define precisely our problem. Following [15| we consider histories 
such that for a given initial configuration a certain integrated current J has flown until some 
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large time T. We then ask how the system behaves during a time interval [t l5 1 2 ] such that 
both t\ and £2 are far from the initial time and the final time T. As shown in [15J this 
is tantamount to computing an effective stochastic process that has j = J/T as its typical 
current and computing the stationary distribution of that effective process. This yields the 
quasi stationary distribution of the original ASEP conditioned on sustaining an atypical 
average flux j. Before proceeding with the discussion of this distribution we make two short 
digressions, the significance of which will become clear immediately further below. 

First we point out that generally the generator of a stochastic process can be written 

H = A-D (4) 

where matrix elements Ac\c — —wc',c < of the off-diagonal matrix A = W+ + W- is given 
by the negative transition rates Wc,c > from a configuration C to another configuration 
C, while the diagonal matrix D = Wq is given by the negative sum of rates out of a given 
configuration, i.e., D C) c — ~J2c ,w c,c- This property ensures that all diagonal elements 
of H are positive real numbers, while all off-diagonal elements are negative real numbers. 
Conservation of probability is encoded in the fact that by construction the matrix elements 
in each column of H sum up to 0. From this property it follows that there is a left eigenvector 
(s I with eigenvalue where all components are equal to 1. This property along with the 
positivity (negativity) and reality condition on the matrix elements guarantee that H is the 
generator of some Markovian stochastic dynamics. 



Generalizing an observation made in 15] we now consider an irreducible matrix M with 
non-negative diagonal and non-positive off-diagonal part. By the Perron- Frobenius theorem 
its lowest eigenvalue, denoted //, below, is real and nondegenerate, and the respective left and 
right eigenvectors (/i|, |/i) have all positive real components /i^ = (C|/i), /i£. = (fi\C). 
Let us define a diagonal matrix A with the components fi^ on the diagonal and consider 
the transformed and shifted matrix 

M' = AM A' 1 - 11. (5) 

It is easy to see that M' has an eigenvalue with a right eigenvector | fjf ) — A| fj, ) and left 
eigenvector (s \ = (/i \A~ l with constant components equal to one. Moreover, by construc- 
tion all off-diagonal elements of M' are non-positive real numbers, while the diagonal of M 
remains unchanged up the constant real-valued shift /1. Hence the matrix M' is the genera- 
tor of a Markov process. The construction ()5]) thus provides a general recipe how to obtain 



a stochastic generator from a rather general family of non-degenerate matrices with non- 
positive real off-diagonal and non- negative real diagonal part. The stationary distribution 
can be written in terms of the lowest eigenvectors of M as 

P* c = nfrg/Z (6) 

where Z = J^ c fic^c = ( A 4 1 f 1 ) * s the normalization that ensures that the sum of all proba- 
bilities is equal to 1. We remark that for symmetric M one has /j,q = ^q. Then M is also 
hermitian and therefore can be seen as the Hamiltonian of some quantum system. We find it 
intriguing that in such a case the stationary distribution of the associated classic stochastic 
dynamics defined by the Markov generator M' is given by the quantum mechanical, i.e., 
quadratic probability distribution of the quantum system. The classical interpretation of 
the presence of both the left and right eigenvectors in this formula is the following : con- 
ditioning on an atypical current up to a final time T requires to produce this current both 



between and t and between t and T — t. Computations made in [15[ show that the left 
eigenvector take into account the first part [0, t] whereas the right one takes into account 
the second part [t, T — £]; there are both a forward-in-time and a backward-in-time effect of 
the global conditioning. 

Next we comment specifically on the current distribution P(J) = Prob[J(t) = J] in the 
ASEP. We introduce the time-averaged current j = J/t, sometimes simply called current, as 
opposed to the integrated current J. Asymptotically the current distribution has the large 
deviation form 

P(J) oc e tf ® (7) 

where f(j) is the large deviation function. Introducing the auxiliary parameter s conjugate 
to the integrated current one obtains the scaled asymptotic cumulant function 

ln(e sJ W) ,. lnV 7 P(J)e sJ W 
fi(s) = lim — i = lim ^ J v ; (8) 

t— >00 t t— >OD t 

which is given by the largest eigenvalue of the modified rate matrix [7|, [l6| 

W(s) = -(e s W + + e~ s W. + W ). (9) 

Here W + = —p J2k s t s k+i anc ^ W- = ~1 ^2k s k s t+i are ^^ e off-diagonal part of the generator 
of the ASEP, corresponding to moves to the right and left resp., and W = J2 k [pn k (l—rik+i) + 
q(l — n^rik+i] is the diagonal part. Thus one obtains the modified rate matrix from the 
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original generator H ([3]) simply by multiplying all off-diagonal elements with e ±s . Notice 
the sign convention that we adopt here to remain consistent with [151 ]: For s = we have 
W(0) = -H. 

The two large deviation functions f(j) and /i(s) are related by a Legendre transformation 

(i(s)=max[f(j) + sj]. (10) 

j 

This construction is in complete analogy to a change between ensembles in equilibrium 
thermodynamics where J would be regarded as extensive variable (here: extensive in time) 
whereas s would be the conjugate intensive quantity. On the one hand one might want to 
study for a given time t space-time histories with a fixed value of J = j CO ndt- Alternatively, 
one could study the transformed ensemble where J is fluctuating, and the conjugate quantity 
s is chosen fixed such that the current has the same mean value 

By construction this "conditional current" is a monotonically increasing function of s and it 
coincides at the point s = with the stationary current of the ASEP. So positive (negative) 
values of s correspond to an atypical current enhanced (reduced ) with respect to the typical 
stationary current of the ASEP. Equivalently one may write 

s(j) = -jJU) (12) 

to obtain the inverse relation between s and j. 

With these notions we can address the question of how the ASEP, conditioned to produce 
an atypical current j during a long period of time, behaves dynamically. We consider this 
question in the ensemble of fixed s rather than fixed j where for a given j the conjugate value 
of s is given by (I12p . As shown in [15J the general construction presented in our digression 
can be used to solve this problem. The rates of the effective Markov process are expressed 
in terms of the left eigenvector of the modified rate matrix (Q corresponding to its largest 
eigenvalue as 

Wg c = Wish J^^y for C^C (13) 

Here Wcc > ^ac are ra t es of change C — V C in the original and in the effective stochastic 
process. The unnormalized stationary state vector of the effective process is given by ([6]) 



and can be written in bra-ket notation as 

P:L(C) = (Ks)\C)(CMs)}. (14) 



The stationary flux (fTTj) of the effective stochastic process is equal to the atypical flux, on 
which the original stochastic process is conditioned. 

From a probabilistic point of view, the transformation from the initial Markov process to 
the conditioned one with the modified matrix M' correspond to a change of measure. The 
expectation value under the new process of any observable A t at time t can be expressed as: 

The existence of the limit and its relation to the conditioning under the current are deduced 
from the large deviation principle satisfied by the current J(T). Such a change of measure 
makes explicit the conjugation relation between the parameter s and the current J(t), similar 
to equilibrium statistical mechanics. 

So, if we succeed to solve the largest eigenvalue problem for the modified rate matrix 
([9]), the effective stochastic process and its respective stationary state can be constructed. 
The aim of this paper is to construct such an effective stochastic process for the ASEP on 
a ring, conditioned to produce a high particle current. The stationary distribution of the 
unconditioned ASEP does not have long-range correlations and the process itself has only has 
hard-core on-site interaction. Conditioning the process on an atypical high flux, one expects 
that the respective space-time particle trajectories (for long periods when this high flux can 
be observed), avoid forming clusters, thus generating an effective repulsive interaction. In 
this paper we find the analytical form of such interaction for very high fluxes. The ASEP is 
chosen for simplicity of presentation, but it will transpire that our results are generalizable 
to non-conservative simple exclusion processes with nearest neighbour hopping, i.e. to a 
class of react ion- diffusion processes. 

We consider first the special cases of two and three particles on a ring of L sites, and then 
generalize to the case of an arbitrary number of particles N. For calculations we use the 
fact that the modified rate matrix for the ASEP belongs for all s to the class of integrable 
models whose eigenf unctions can be found using the Bethe Ansatz. 

We stress that three different connections of classical stochastic processes to quantum 
mechanics are described in this work. They are independent from each other and have their 
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own meaning. The first one presented in the previous section is the quantum Hamiltonian 
formalism to describe the evolution of the probability distribution. It arises from the Markov 
property of the stochastic process, in particular from the linearity of the temporal evolution 
of the process. However, the classical probabilities are linear in the right eigenvector, while 
quantum mechanical expectations are quadratic. The second connection, which is new, stems 
from the effective dynamics where the conditioned expectation values are now quadratic in 
the eigenvectors (see (I14p ) as in quantum expectations. This fact is independent of the 
previous quantum formalism. Finally, the determinantal free fermion structure described 
in the next sections for a large current is also a new feature of the present regime of the 
exclusion process and is not a direct consequence of the two previous facts. 

III. TWO PARTICLES ON A RING 

The eigenfunction of the modified rate matrix W is found using the coordinate Bethe 
ansatz (for the present context see e.g. [9l. Il6l.l22j]) 



\f)= Yl (zi 1 z?+A 21 z?z^)\n 1 n< 2 )= ^ Y nin2 \ ni n 2 ) (16) 

l<nl<n 2 <L l<nl<n 2 <L 

This eigenvector has eigenvalues A = pe s /zi + pe s / z 2 + qe~ s Zi + qe~ s z 2 — 2 and the ampli- 
tude A 21 is given by requirement that Y nin2 is an eigenfunction for particles being nearest 

neighbours 

A 21 = _ ^ + g e-^ 2 -, 2 

pe s + ge s z±z 2 — z\ 

Moreover, periodic boundary conditions require 

z 2 = A 21 . (18) 

The translation invariance of the problem together with the constant sign of the compo- 
nents of the ground state yields Y nin2 = Y ni+1>n2+1 and therefore z\z 2 = 1 for the ground 
state. It is convenient to parametrize z 2 = e n , z\ = e~ n . Then (fTB"|) is a transcendental 
equation which has L solutions for 7. The amplitude corresponding to a configuration with 
two particles at distance / follows as 

Y n , n+l = 4 + A 21 z[ = 4 + 44 = J* + e^ L ~ l \ (19) 
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We have to determine that value of 7 that gives the largest eigenvalue. From the fact that 
the largest eigenvalue A = (pe s + ge _s )(e* 7 + e -ry ) — 2 = fi(s) is real, it follows that 7 is either 
real or imaginary. 7 as function of s is determined by solving ( TT71) with substitution (TT81) . 
For s close to 0, it turns out that real 7 corresponds to s > and imaginary 7 corresponds 
to s < 0. For large e s ^> 1, 7 is real and the Bethe amplitude Y n>n+ i ( fT9l) of the largest right 
and left eigenvectors is given by 

F n , n+/ = 2sin7r//L + 0(e- s ) (20) 

The origin of the above formula will be explained in the section |Vl 

IV. THREE PARTICLES 

The case of two particles was treated as an introductory pedagogical example. In this 
section we study the case of three particles in more detail. 

A. Bethe ansatz 

The eigenfunction of the modified rate matrix W for A^ = 3 particles on a ring is given 
by the Bethe ansatz 

W)= Yl ^2 A ° z T(i) z 7(2) z T(3)\ n i n 2 n 3) = Yl Yn in2 n 3 \nin 2 n 3 ) (21) 

l<nl<n2<n-j<L cr l<nl<n2<ri3<L 

where a are the 3! = 6 permutations of indices ni, n 2 , nz- This eigenfunction has eigenvalues 

3 3 

A = pe s J2 l /*i + <? e ~* J> - 3 (22) 

and the amplitudes A a are given by requirement that (p is an eigenfunction for particles 



being nearest neighbours [9|, lla, |22[). This yields e.g. 

(23) 



Ajik _ aji _ pe s + qe s ZiZj - Zj 



Aijk otij pe s + qe~ s ZiZj - z t 

and similar relations for any other permutation of two indices. From the periodic boundary 
conditions it follows that 

f[ (-1)2* = zl = 17 peS + qe ~ Sz * Zk - Zk (24) 

J--L a ^ {}pe s + qe- s ZiZ k -Zi 
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for all k. 

As in the previous case N = 2, from translation invariance for the ground state compo- 
nents we have Y nin2n3 = Y ni+ i tn2+ i tna+ i we get z\z 2 z 3 = 1, so we have only two unknown 
parameters. They must be found from the Bethe Ansatz equations ( 124|) and the requirement 
of maximization of the eigenvalue A. Using ( 1231) . ( EH1) . and z\ z 2 z 3 = 1, the Bethe amplitudes 
Y nin2n3 for 7i\ = n,n 2 = n + l,n 3 = n + I + d may be brought to the form 

Y n ,n+l,n+l+d = A 123 [z 2 Z 3 + Z 3 Z ± + Z X Z 2 Z 3 } + A 2i3 [1 -H- 2] (25) 

where A 213 = — a 2 i/«i2 = —(pe s + qe~ s ziz 2 — z 2 )/(pe s + qe~ s z 1 z 2 — Zi), and [1 <H> 2] stands 
for the substitution z\ -H- z 2 in the expression inside the first square brackets in (125]) . 

B. The limit of large conditioned current 

Consider the limit e s ^> 1, which corresponds to a dynamic restriction on an atypically 
large current. Treating e~ s as a small parameter, we can write the Bethe equations in the 
form 

.L 



~k 



1 + e- V 1 I Y. z i ~ 3zk j + °( e ~ 2 *) ( 26 ) 

We look for the solution for Zk in the form Zk = e %lk . We choose 7^ as 

lk L = 2(k - 2)it + a k e~ s + 0(e~ 2s ). (27) 

This choice leads to the maximal eigenvalue we are searching (basically, we have to pick up 
the set of z k with the largest real part). Note that in the following we shall also assume the 
large L limit L> 1. Solving self-consistently the Bethe equations, we obtain after some 
algebra a k = S'jk/p, so that 

2(fc-2W 3 __.\ , „,_„,, 



■*=^-L^\?-VL e ~) +0(e '' <28) 

The respective maximal eigenvalue of the modified rate matrix is given by 

3 3 

H{s) = pe s J2 z k l + ^~ S J2 Zk ~ 3 = P eSJ ° + °( 1 ) + °( e "') ^ 29 ) 

fc=i fc=i 

where 



Jo = J2cos 2 -^^ = 3-% + 0(l/n (30) 



fe=i 
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The respective magnitude of the enhanced current is dfi(s)/ds = pe s J . So indeed the choice 

e s 3> 1 corresponds to large atypical current. 

Substituting the solution of the Bethe equations ( 1281) into (1251) . we obtain the right eigen- 

function F123, which corresponds to a configuration with 3 particles, separated by distances 

h, hi h where Yl lj = L. The result is (up to an overall constant) 

3 3 

^i23 = 2^sin 7 4-e^sin7^cos7/ fc + 0(e- 2s ), (31) 

fe=i k=i 

where 7 = ^ (1 — 3e~ s /L), or, alternatively, as 

Y 123 = 8 sin gZi) sin (|z 2 ) sin (^/ 3 ) + 0(e~ s ). (32) 

The left eigenfunction turns out to be the same as the right eigenfunction, which can be 
demonstrated either from the Gallavotti-Cohen symmetry or by the following argument. The 
transposition of the matrix ( W(s) ) results in the exchanging pe s •<=>■ qe~ s in the original 
matrix W(s). Denote the Bethe roots of the transposed matrix as Tk- In can be verified 
that the Bethe Ansatz equations and eigenvalue of the transposed matrix coincide with the 
old ones (123|) . ( 122]) . after substitution r k — >■ l/z k for all k. Under these substitutions, the 
amplitudes of the left eigenfunction become Y^ FT = Y i A a z~J 1 j lz~J^ j z~J^y However, note 

IT 

from (1281) that z\ = z^ 1 ,and z 2 = 1, up to corrections 0((e _s ) 2 ). Therefore, by relabelling 
the roots Zl O z 3 we obtain Y™ FT = ^ A^f^fah = Y nl GHT - 

(7 

With that in mind, we are in a position to construct the effective stochastic dynamics, 
defined by (!T3|) conditioned on large flux. 

To this end, notice first that the stationary probabilities to find the configuration of 
particles at distances h,h,h, is given by (JT4J) 



P cS (h,l2,h) = (Y 123 ) 2 /Z (33) 

where Z = ( J^ fi l2 Zjj (F123) ) = 3L 2 is the normalization factor. With ( !32l it is readily 
checked that 

• (A) the configuration with h = h = h = L/3 is the most probable one 

• (B) all configurations with any single finite distance /&. i.e., where h/L ^C 1, and 
arbitrary other distances lj have vanishing probability 
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Indeed, this is consistent with the intuition that the configurations where two particles are 
coming near each other (close on lattice scale) should be avoided in the effective stochastic 
dynamics, because it would tend to cause a current reduction due to hard-core repulsion. 

The rates for the effective stochastic dynamics are given by ( TT3l) . Namely, the rate for 
hopping to the right of the second particle at position X2 in the configuration with consecutive 
distances fa = x i+ i — Xj, is given by 



W, 



cir 



h+l,h-l,h\hhl3 ~ P e 



sin (f (fa + 1)) sin (j (fa - 1)) 
sin (lh) sin (f fa) 



0(e 



(34) 



(remember that Wcc is rate of change C — )■ C"). At short distances fa/Z/ < 1 or fa/ 'L ^C 1 
the major contribution to this expression in square brackets is calculated by an expansion 
sine « e. We obtain that the rates of hopping for a two particles at small distances d , 
d/L <C 1 are given by 

(35) 



-cff 
d , d+l.. 


^ pe s 


d 
_d + l_ 


, d = l,2, 


-cff 
..d+l d 


pa pe s 


'd+l' 

d 


, d = l,2, 



(36) 

At the same time, the rates of hopping to the left in the effective dynamics will have the 
prefactor qe~ s , and in the limit e~ s <C 1 can be neglected. We see that the rates (1361) 
driving two particles apart are larger than those driving them together (135]) . and therefore 
the configurations with two particles forming a cluster are difficult to reach. However, in 
any finite system L < oo, even in the limit e~ s — > 0, configurations with particle distances 
(/j = 1) are possible. 

It is instructive at this point to compare the probability of forming of a cluster in the 
model with the effective dynamics and in the original ASEP. In the model with the effective 
dynamics the single-cluster probability is 

2 8 7T 6 



pcff 



OI l 1 =l,l 2 =l,l z =L-2 



/(3Z/ 



^ 



L s 



for L > 1, 



(37) 



and the probability of a cluster of two particles is 



L-3 



x-ff 



3L 2 



E^i 



h=2 



W3=L-i-^2 6 ^L(sin 4 x)/(L 2 



247T 2 

~1F 



for L > 1. 



(38) 



On the other hand in the original unconstrained ASEP with 3 particles on a ring, the 
probability to form a single cluster is P^ SEP = £/( 3 ) — > Q/L 2 for L>1, and the probability 
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to form a cluster of two particles is P 2 ASEP = L(L — 4)/( 3 ) — > 6/L for L>1. So we see that 
the clustering is strongly suppressed in the effective dynamics, with respect to the unforced 
ASEP, which is natural to expect in a system conditioned to produce a large current. On 
the other hand, the stationary probability of equidistant state h = h = h = L/3 (the most 
probable state in the effective dynamics) 

P C °U = Y? i=h = h =L/J^) = 9/L 2 (39) 

is only 4.5 times larger than in the unconstrained ASEP P^Ef d p = f/( 3 ) — > 2/L 2 for L ^> 1. 

Finally, we can define an effective potential U, felt by the particles evolving according to 

the effective stochastic dynamics, in obvious analogy with equilibrium thermodynamics, by 

(3U{C) = -lnP cff (C) = -2ln(^i(s)|C), (40) 

where (3 is the effective inverse temperature. Consider first the case of two particles, where 
the amplitude (/xi(s)|C) for the configuration of two particles at distance I is given by 
Y 12 = simrl/L + 0(e' s ). Hence 

71-7 

PU(1) = -2 In Y 12 = -2 In sin — + const (41) 

shown on Figfl] (solid line). By numerical evaluation of the Bethe ansatz equation we have 
also determined the effective potential for finite s. Analogously the effective potential for 
three particles takes form 

(3U = -2\nY 123 (42) 

where 1x23 is given by formula ( 13T|) . Note that below we shall use the freedom in defining 
the potential by eliminating a constant, connected to the normalization of the P eS (C). 

C. First order corrections in e s to the effective dynamics 

Using (1311) . one can obtain first order correction to the effective dynamics. E.g., for the 
hopping rates to the right ( |35l) we get, assuming d/L <^i 1 

rcS , x „ d {_ e~ s l + 5cos27ru 



W'^Js) = pe*— (1 - d{d+l)1 _ m2m + O(e-) ) , (43) 

where u = h/L, for a configuration with two particles at small distance d. 
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FIG. 1: Effective Potential U felt by two particles on a ring of L = 100 sites. The short-distance 
configurations (l/L = 0, 1) are unlikely. The solid, dashed and dotted lines correspond to e~ s = 
0, 0.5, 0.9 respectively. 



For generic configurations with particles separated by distances which scale with L, we 
get, by expanding the generic expression ( 1T3l) with substitution ( l3TT) in Taylor series with 
respect to l/L and e~ s and keeping the first order terms 



W h+l,l2-l,h\hl2h^ 



pe 



27 r (cos^-cos^) 2tt $(£)-$ (£) 



27dfc 

L 



^ 2 ELi* 



(44) 



where $ (w) = ^ [(1 -+- 3m) sin27ru]. The last term in the brackets gives the contribution due 
to finite s. The relative effective potential differences is given the ratio between the middle 
and the last term in the brackets, 



AU(s) - AU{s -► oo) _ e" s $ (|) - $ (£) 



L 



cos 



2?rh 



COS 



27T« 2 ' 



(45) 



AJ7(s -> oo) 

As the system size L becomes infinite, the first order corrections to the hopping rates due 
to finite s survive for configurations with finite (on lattice scale) distances as (1431) shows, 
because it does not have l/L dependence in it. On the contrary, for configurations with all 
inter-particle distances scaling with system size < U/L < 1, for L — > oo, the corrections 
due to finite s disappear, see (1441) . (145]) . Another example of this can be seen in FigJU 



17 



V. ARBITRARY NUMBER OF PARTICLES IN THE LARGE CURRENT 

LIMIT 

A. Largest eigenvalue and associated eigenvector 

The arguments, discussed in previous sections, are readily generalized to the case of an 
arbitrary number of particles N. We shall treat the large current limit e s — > oo only Note 
that unless stated otherwise, the system size L is arbitrary. Note also that in the limit 
e s — > oo , to avoid the divergences, the modified rate matrix ([9]) has to be rescaled by the 
factor e~ s as e~ s W(s). 

For finite s the eigenfunction of the modified rate matrix is given by the Bethe Ansatz 
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22|] 



<p)= Yl Yl A ° z ?(i) z *(2)--- z a(N)\ n i n 2--- n N) = Yl Yi2..N\nin 2 ...n N ) 

I<?il<n2<n3<l/ crSS'jv l<nl<n2<n3<L 



(46) 



(47) 



where the coefficients A a are given by 

A...g(jj)... _ _Oji L _ pe s + qe~ s ZjZj - Zj 
A..ij.... ocij pe s + qe- s ZiZj - z t 

and the Bethe roots Zj satisfy 

T7 (-1) °* = z L k = (-l)^ 1 f] P&S + qe ~ SZlZk ~ Zk (48) 

. ,.,, ot ik ±J - pe s + qe s ZiZ k - Zi 

for arbitrary N. 

In the limit e s — > oo the solution of the Bethe equations corresponding to the largest 
eigenvalue of the modified rate matrix has the form Zk = exp(z7fc), where 

lk = ^-A k = l,2,...N (49) 

generalizing (128]) . The maximal eigenvalue of the rescaled modified rate matrix e~ s W(s) 
and the respective rescaled current are given by 

N 

lim e~ s n(s) = p^ zV 1 = J (50) 

fc=l 

where 

^o=P> cos =P^— 7 -■ (51) 

*--' L sin t 

fc=l L 
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The coefficients A a are simplified to 



A a = sgn(a) 



(52) 



where sgn(a) denotes a signature of permutation a . The limiting expression for Y\2..n can 
then be rewritten in the equivalent form 



Y 



12. .AM 



^2 sgn(a)zl 



°"(™i)^°"( n 2) a{n N ) 



Z ...Z 



N 



(53) 



ctGSjv 



In the following we drop the subscript indicating that we are taking the limit of large 
s. Using the Leibniz formula for the determinant, we can rewrite the above expression as a 
Slater determinant 



Y 



12..N 



det 



,.»i 



-2 



zT 



y "2 

'2 



™]V 



y n N 



z N z N ... z N 



(54) 



For the choice of Bethe roots (I49D this determinant takes the form 



Y 



12..N = exp 



A' 



{N-l)"£i 



ITTXj 



3=1 



n (■ 

l,k 
l<Kk<N 



2i-KXi. 2inx 

! " — e 



") 



(55) 



where x k = n^/L, which becomes simply 

Y 12 ...N = 2(a) J] 



sm7r- 



n k -rii 



(56) 



l,k 
l<Kk<N 



The amplitudes Y\2...n at this stage do not have a probabilistic meaning and are not nor- 
malised. 

For N = 2,3 we recover the expressions (I2"0j) . (13"2"j) up to finite s corrections. By ex- 
tending the arguments given after (132]) we readily verify that the left maximal eigenvector 
is given by the same expression fl56|) . Consequently, the stationary probabilities of particle 
configurations in the effective dynamics are equal to the normalized square of the amplitudes 
PJ, 



12. ..N 



K 



h n 



.2 n k-n L 
sin it . 



l,k 
l<Kk<N 



where 



K N T , = 2- N ^ N -^L N 



(57) 



(58) 
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is the free fermion normalization factor [18(. Notice that in this result, which is exact for any 
system size L > 1 the particle number N is subject only to the exclusion rule < N < L. As 
detailed in section |VU] further below, this formula is similar to the distribution of eigenvalues 
in the circular unitary ensemble of random matrices. 

Let us analyze our main results ( 1561) . ( 1571) . The stationary probabilities depend on the 
scaled macroscopic distances x k . For N fixed and finite they satisfy the following properties 
in the thermodynamic limit L — > oo: 

(A) max Y12...N is reached for \xk+i — Xk\ = — for all k (59) 

N 

(B) Y 12 ... N = for 11 (x k+1 - x k ) = (60) 

fc=i 

The property (A) indicates that the most probable configuration of the system is an equidis- 
tant configuration. The property (B) is evident from (15 6p and implies that particles cannot 
be found with finite probability at distances which are finite on lattice scale (or in fact in 
any non-macroscopic distance o(L)). 

B. Effective stochastic process 

Now we describe the effective stochastic process corresponding to the limit e s — > 00. The 
effective hopping rates are given by ( TT3|) : The move of the fc-th particle, located at position 
nk = Lxki to consecutive position n k + 1, provided it is vacant, has the rate 

hm e-V 1 ^ = ^ = II ^^V- ( 61 ) 

s->oo Y12..N ,Jr sm7r k T '■ 

where the product is over I and k as above. The stationary state probabilities with respect 
to the above rates are given by (|57|) . 

The effective potential of interaction, given by ( l40l . has the form of a long-range potential 

U(xi,x 2 ,.-,x N ) = -2 ^J log I sin ir(xi - x k )\. (62) 

l,k 

Note that the summation in ( )62|) extends over all % ^ j. Note also that due to the ring 
geometry and the fact that | sin(7r ± nc)\ = \ sin7rc|, the differences x\ — x k in f )62|l can be 
substituted by the smallest distances between particles I, k on a ring. As for the stationary 
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distribution of the particles ( 157)) . this effective potential U is similar to the potential that 
appears for eigenvalues in random unitary matrices that perform Dyson's Brownian motion 
over the unitary group U(N) (see details below in section IV CI) . 

For any finite number of particles N/L <C 1, the long range effective potential U attains 
in the thermodynamic limit L — > oo singularities for any configuration with finite distances 
(on lattice scale) between the neighboring particles. For the corresponding hopping rates of 
the k-th particle we find to leading order in 1/L 

lim e-y 1 ^ = 1 + y^- lnY 12 .. N = 

= 1 + j ^cot7r(x fc -xi). (63) 

This is obtained by expanding ( IBT)) in a Taylor series with respect to 1/L. Such dynamics 
correspond to an enhancement of the hopping bias whixh vanishes at macroscopic distances, 
but which is strong at microscopic distances. 

It is instructive to compare some configurational probabilities in the ASEP with the 
effective model for finite fraction of particles. For the half-filled lattice N/L = 1/2, the 
probability of the equidistant particle state 0, 1, 0, 1, ... in the ASEP (without conditioning) 
is given by P^sep ("0101..) = 2/( L / 2 ) ~ 2" L for large L. In the effective dynamics, the 
probability of the same state (the state of largest probability ) is expressed, after some 
algebra, in terms of the normalization factors ( 1551) as P eff (..0101..) = 2Kl/2,l/2/ 'Kl/2,l — 
2 1_ "2 . So we see that 



P eff (..010L.) = ^/P^p S *("0101"), (64) 

up to terms which grow in size slower than exponentially. 

This calculation can be easily extended for arbitrary finite filling fraction of the form 
N/L = vp = 1/2,1/3,... The ASEP weight of the equidistant configuration for L 3> 1 
is P^p st (100...100...100..) = ^(zij -1 ~ e^ 111 ^^ 1 -^ 1 ^ 1 -^, and the respective 

l/v F 

weight under the effective dynamics is given by P cS (uf) = Kl Uf ,lv f /{vfKlv f ,l) = Vp F ~ ~ 
Q Lv F \nu F C om paring the two probabilities, we obtain, assuming L 3> 1: 

p%™f\ioo... 400^1 00..) 

1/VF - p L(l-v F )ln(l-u F ) / fi KN 

Moreover, using the connection to the quantum free fermion model at zero temperature 
(see Conclusions), various correlation functions for the effective stochastic dynamics can 
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be obtained from the well-known free fermion quantum expectation values. The particle- 
particle correlation function in the effective dynamics {fikfik+m} is equal to the ground state 
expectation value ((1 — crf)(l — cr f +m ))/4 in the free fermion problem and is given in the 
thermodynamic limit L — y oo, N/L — y p by 

«ww>- = p 2 - A 5 ^ («) 



see e.g. (17|| . Note that in the unconstrained ASEP there are no correlations in the stationary 
state (h k h k+m ) asep = (nk)(h k+m ) = p 2 for any m > 0. For half-filling p = 1/2 and nearest- 
neighbour correlations we obtain (ll) cff = (00) cff = 1/4 - 1/tt 2 « 0.149, and (10) off = 
(01} cff = 1/4+1/tt 2 « 0.351 for the effective dynamics while (11) = (00) = (10) = (01) = 1/4 
for the original ASEP. This again shows that the effective dynamics of ASEP conditioned on 
high current avoids building clusters, as is expected. Interestingly, for p = 1/2 the effective 
particle particle correlations fl66|) at even distances m = 2n coincide with those for ASEP, 
(nfc?7,fc + 2„,) cff = (nkfik+2n) ASEP = 1/4 while the odd-distance correlations decay algebraically 
(n k n k+{2n+1) ) cS = 1/4 - l/(7r 2 (2n + l) 2 ). 

C. Circular unitary ensemble and Dyson's Brownian motion 



Exclusion processes on Z are already known to share common features [Hi, ll9|, [20|] with 
random matrices, such as the determinantal structure, the Airy fluctuations and the Tracy 
Widom distribution for different quantities in the TASEP. Howver, much less is known about 
the connection between the ASEP on a ring with periodic boundary conditions and random 
matrix ensembles. 

The circular unitary ensemble is defined as the Haar measure on the unitary group U{N). 
This group is defined as the set of complex matrices U such that XPU = I. This Lie 
group can be equipped with a uniform measure, called the Haar measure. Almost every 
matrix U G U(N) can be diagonalized and the Haar measure induces a new measure on the 
eigenvalues, which contains a determinantal-like interaction between them. All eigenvalues 
of a unitary matrix lie on the unit circle U and can be written as e l9k where < 9 k < 2tt. 
The density probability of eigenvalues is 

PcaB{0 1 ,...,0 N ) = j- J] |e" fc -e"'| 2 =;Jr II sin2 (^) > ( 67 ) 

N Kk<KN N Kk<KN ^ ' 
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which is precisely the formula (157"]) . up to the fact that the n,j/L are rational numbers due 
to the discrete structure of the underlying lattice. 



This analogy can be developed further if one considers Dyson's Brownian motion 2l| on 
the unitary group U(N). The infinitesimal evolution of the process Ut on U{N) starting at 
the identity matrix is given by (Ito form): 

1 



dU t = idH t Ut - -U t dt 

Li 



and can also be seen as the exponential 



U 



JH t 



(68) 



(69) 



where H t is a random hermitian matrix whose independent entries perform independent 
Brownian motions. The independent diffusions of the coefficients induce interacting diffu- 
sions on the eigenvalues e l9k ^ of the matrix U t : 



d0 fc = £cot(^M 



(70) 



»(*) 



where the B\ are N independent Brownian motions. Up to a global common drift term, 
the interaction term is similar to the one obtained in (163]) for the ASEP conditioned to 
produce a large current and it derives from the same long-range potential U as defined in 
(I62p . The invariant measure of the Dyson Brownian motion over U(N) is precisely the Haar 
measure ( 167]) . 

Although the similarity is striking, a deep understanding of the connection between both 
models and the existence of further similarities (dynamical fluctuations, corrections in e~ s ) 
are still missing. 



D. Vandermonde determinant limit and clustering probabilities 



Let us consider a special case of finite number of particles N, which form a single cluster, 
where all distances between particles are small, i.e. \xu — xi\ <^C 1 for all k, I. In this case the 
amplitude of the wave function (156]) becomes proportional to the Vandermonde determinant 



Y 12 ...N OC (Xk 



Xl 



det 



l<l<k<N 



1 X\ 
1 X 2 

1 X N 



X 



JV-1 

1 



x 2 



X 



JV-1 

JV 



(71) 
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The respective stationary probabilities of the configurations are proportional the square of 
the Vandermonde determinant, 

P xl x 2 ...x N oc (x k - x{f (72) 

l<2<fc<7V 

For instance, the ratio of probabilities of two similar single cluster configurations C, C where 
C differs from C by rescaling all particle-particle distances by a common factor r, namely, 

\x' k — x[\ = r\xk — xi\ for all k, /, is given by 

peff 

"V =r iV(iV-i) ) (73) 



peff 

1 XlX 2 ...X N 



which grows much faster than exponentially with number of particles N (an exponential 
growth would be expected for uncorrelated particles, restricted to smaller volume). E.g. 
single cluster configurations of 3 particles (10 particles) are 2 6 times (2 90 times) less likely 
than respective similar configurations of double size r = 2. 

E. Gap in the spectrum of the effective Hamiltonian 



The expression (1561) gives the components of the maximal eigenvector of the effective 
Markov rate matrix W cS which corresponds to the stationary state. Let us find the next 
largest eigenvalue of W which determines the longest relaxation time. Following our 
introductory exposition (see also [15|) all eigenvalues of W eS are given by Aj(s) — fi(s) 



where Aj-(s) are eigenvalues of the modified rate matrix W(s) (Q. In the present case the 
next largest eigenvalue of W(s) is given by the set of Bethe roots for the largest one ( )49l ). 
where the last root zjy = exp (i^j^-) is replaced by the root z' N = exp (i^^ )■ Such a choice 
corresponds to a minimal modification in the set (H9l) and it is inherited from the underlying 
free fermion quantum system. 

Computing the difference between the eigenvalues A = /x(s) and A' given by the set (I49p 
and the new set, we obtain to leading order in e s the real part of the spectral gap of W , 
as 

Re lfi(s) — A'] ~ pe s (cos7r(p — — ) — cos7r(p + — )) 

LJ LJ 

TX 

= 2pe s sin — sin up (74) 

Lj 
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where p = N/L is the particle density. Since our choice of Bethe roots corresponds to 
e s — > oo limit, we see that for large L the gap in the rescaled effective Hamiltonian is 
inversely proportional to the system size e~ s W cG = 0(1/ Li). Note that in the unconstrained 
ASEP the real part of the spectrum gap scales as 0(1/L 3 / 2 ). We conclude that the long- 
range interaction of the effective process leads to a marked reduction of the longest relaxation 
time (which is the inverse spectral gap). 

VI. CONCLUSIONS 

We have shown that conditioning the ASEP on a long period of strongly enhanced current 
can be described in terms of an effective exclusion process with long range interactions. 
Numerical evaluation of Bethe ansatz equations and perturbation theory in e _s show that 
the explicit asymptotic results that we have obtained for very large fluxes are rather robust 
with respect to the value of s, at least for finite number of particles N/L <C 1 and large 
system sizes L. Indeed, as can be seen in Figfj], the effective two-particles potential U for 
e _s < 0.5 is difficult to distinguish from the limit e _s = for L = 100. For L = 1000, the 
differences from the limiting case would become visible in FigJT]only at values e _s > 0.95 
and higher (curves not shown). The case e _s = 1 (s — 0) corresponds to the usual ASEP 
dynamics on a ring without conditioning and presents a singular limit. 

The hopping rates corresponding to the limiting effective dynamics have been computed 
explicitly, see (l6"Tj) . They can be interpreted as arising from an equilibrium potential that 
depends only on the interparticle distance. We have computed various space-time properties 
of the dynamics conditioned to carry a large current for a long time, i.e., of the quasi- 
stationary distribution. Density correlations decay with a power law and have a periodic 
modulation 0661) . Temporal correlations decay with a dynamical exponent z — 1 rather than 
with the KPZ exponent z = 3/2 that characterizes the typical behaviour of the ASEP. Hence 
in the extreme-current regime the process is in a different universality class, indicating a 
phase transition as one enters into the regime of strongly non-typical behaviour. Non- 
analytic behaviour of physical observables in the regime of strong non-typical behaviour has 
been shown to occur also for atypical small currents and in the zero-range process for high 
currents. 

Our approach highlights a surprising analogy between the classical stochastic problem 
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we have started from with an intrinsically quantum problem of N free fermions on a periodic 
lattice. This analogy goes far beyond the quantum formalism used in section |TT] to describe 
the stochastic classical dynamics. In fact, the expression ( I54p . and, correspondingly, also 
( |56l) is a Slater determinant describing the wave function of a system of free fermions, the 
so-called XXO model, 

L 

HxXO = ~ Yl { a n<+l + ^n+l^n ) ( 75 ) 

n=l 

with periodic boundary conditions L + 1 = 1. The probabilities of the classical particle 
configurations under effective dynamics given by P cG (C) = (/ii(s)|C)(C|/xi(s)), coincide, up 
to terms of order e~ s , with quantum expectation values for the free fermion problem because 
of the presence of both left and right eigenvectors, whereas the unconditioned probability 
involves only the right eigenvector. The reason why the free fermion Hamiltonian has a 
common eigenvector with the modified rate matrix in the limit e _s — > is simple: In this 
limit, the modified rate matrix becomes e~ s W(s) = J2 n =i { a n (T n+i)i which commutes with 
XXO Hamiltonian. 

Our asymptotic results are not limited to the ASEP and they are also not limited to 
conditioning on the current. One would obtain the same spatio-temporal organization of the 
particle system in the limit of large activity, i.e., counting the number of jumps irrespective 
of direction. In that case the modified rate matrix (J9J) would take the form W(s) = e s (W + + 
W-) + Wq- In this case the limit e~ s — > would directly lead to the XXO Hamiltonian (1751) . 

Similarly one can study non-conservative exclusion processes where particles can be cre- 
ated or annihilated, see e.g. the 16-parameter family of reaction diffusion processes defined 
in [5|, |6J. This includes diffusion-limited annihilation, the contact process and other well- 
studied processes. For all these processes the quasi-stationary distribution for large current 
or activity is the same as for the ASEP computed here. Hence our asymptotic results apply 
to a large family of interacting particle systems. 

The connection to the Circular Unitary Ensemble for s — > oo and the already known 
connections of the ASEP with the Gaussian Unitary Ensemble raise the question of further 
connections with random matrices. In particular, it would be interesting to perform the first 
corrections in e~ s of the effective probability P cS , in order to check if it corresponds to the 
corrections of some potential on the unitary group U{N). 
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